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FOREWORD 

This  report  describes  work  directed  toward  estimating  the  effect  of  multiple 
Ejection  Rack  (MER)  flexibility  on  store  ejection  conditions.  This  work  was 
sponsored  by  the  Naval  Air  Systems  Command  under  AIRTASK 
A3200000/009B/3F232000.  The  report  was  reviewed  by  Dr.  Thomas  A.  Clare,  Head, 
Flight  Dynamics  Group,  and  Russell  D.  Cuddy,  Head,  Aeroballistics  Division. 
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Structural  flexibility  has  long  been  suspected  as  being  a  factor 
affecting  the  ejection  phase  dynamic  response  of  stores.  After  hook 
opening,  in-carriage  loads,  associated  with  a  released  store,  and  an 
ejection  recoil  force  deflect  the  MER  beam  away  from  its  initial  position. 
This  in  turn  leads  to  a  lower  ejection  velocity,  lower  ejection  pitch 
rate,  and  an  induced  rolling  moment.  A  theoretical  structural  dynamic 


00  uSTi  1473  COITION  or  «  NOV ••  II  OMOLCTt  UNFI  ASSIPf m 

*/n  ,  *  untmsamED  _  _  u 

|  SECURITY  CbAlllf  ICATION  Of  Tni»  RaOC  f»hdn  line  United) 


jMjTv  Cl*S4lflC*T|'3>;  Or  Twts  P*G£(wh*n  r-**  &>••**) 


model  is  based  upon  "Bernoulli"  pitch  and  yaw  deflections  and  torsionc.l 
rotation  of  the  HER  beam;  wing,  pylon,  hanger,  shoulder  pad,  and  ejecxor 
unit  deformations  are  neglected. 

Detailed  development  for  the  present  report  is  confined  to  the  pitch 
plane  bending  and  store  pitch  dynamics  case.  For  small  pitch  angles 
a  linear  set  of  kinematic  and  dunamic  equations  for  the  store  ejection 
phase  is  obtained.  Sample  computations  are  given  for  two  stores  and 
constant  aircraft  puli-up  rates.  A  4-g  pull-up  produces  the  greatest 
structural  deformations.  However,  the  pull-up  maneuver  alone  affects 
the  ejection  pitch  rate  and  velocity  to  a  greater  extent  than  the 
structural  deformation.  A  break  even  point  occurs  for  smaller  pull-up 
rates.  For  the  pitch  bending  case,  flexibility  probably  accounts  for  no 
more  than  a  105  deviation  from  the  rigid  case.  Flexibioity  affactj  the 
ejection  angular  rate  more  significantly. 
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the  MER  beam:  wing.  p> Ion.  hanger,  shoulder  pad.  and  ejector  unit  deformations  are 
neglected. 

Detailed  development  for  the  present  report  is  confined  to  the  pitch  plane 
bending  and  store  pitch  dynamics  ease  For  small  pitch  angles  a  iinear  set  of 
kinematic  and  dynamic  equations  for  the  store  ejection  phase  is  ootained.  Sample 
computations  are  given  for  two  stores  and  constant  aircraft  pull-up  rates.  A  4-g 
pull-up  produces  the  greatest  structural  deformations.  However,  the  pull-up  maneuver 
alone  affects  the  ejection  pitch  rate  and  velocity  to  a  greater  extent  than  the 
structural  deformation.  A  break  even  point  occurs  for  smaller  pull-up  rates.  For  the 
piUh  bending  case,  flexibility  probably  accounts  for  no  more  than  a  UK?  deviation 
from  the  rigid  case.  Flexibility  affects  tne  ejection  angular  rate  more  significantly. 
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I.  INTRODUCTION 


Stoietura!  flexibility  has  long  been  suspected  as  being  a  factor  affecting  the 
ejection  phase  dynamic  response  of  s*ores  released  from  Multiple  Ejection  Racks 
(MER).  Reference  1.  which  provides  a  good  summary  of  the  problem  areas  of  store 
separation,  contains  several  qualitative  references  to  flexible  MER  effects. 
Reference  2  considers  'he  effect  of  rack  deformation  on  ejection  velocity  from  a 
simple  energy  consideration  for  a  uniform  cantilever  beam  with  a  half  sine  pulse 
ejector  force.  Reference  3.  considers  a  simple  one  degree  of  freedom  model  for  a 
bomb  rack  and  separating  store  with  step  function  ejector  force. 

The  objective  of  the  present  report  is  to  approximate  the  structural  dynamic 
response  of  the  MER  ejector  beam  and  the  ejection  phase  of  store  separation  in 
more  detail.  Detailed  developme  it  will  be  limited  to  an  aircraft  pulling  up  at  a 
constant  rate,  pitch  plane  bending,  and  pitch  plane  motion  of  the  released  store. 
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II.  STRUCTURAL  DYNAMIC  MODEL  OF  /•»  MER  EJECTOR  BEAM 


I'lBUic  I  Ih  (i  schematic  of  a  MLR  7  or  9  ejector  brum.  The  following 
sequence  of  events  describes  most  of  the  pertinent  ejection  phase  phenomena, 
associated  with  a  flexible  rack,  which  occur  after  hook  opening: 

(1)  The  In-eurriage  weight,  acceleration,  and  aerodynamic  loads,  associated  with 
a  released  store,  are  removed  from  an  Initially  deflected  MUR  beam, 

(2)  The  ejection  force  displaces  the  store  until  the  ejector  foot  Is  extended 
about  3.25  Inches,  Simultaneously  a  recoil  force  acts  on  the  ejector  unit, 

(3)  The  MUR  beam  deflects  away  front  Its  initial  position,  This  changes  the 
position  of  the  ejector  foot  relative  to  the  store.  For  a  single  store  or 
long  timvlntcrval  string  release  flexibility  contributes  to  a  lower  ejection 
velocity,  lower  ejection  pitch  rate  and  an  Induced  rolling  moment,  For  a 
ripple  or  salvo  mode  the  beam  Is  Initially  in  motion  and  release 
phenomena  are  very  much  dependent  on  the  release  Interval. 

In  Figure  I  are  shown  coordinate  systems  used,  x,  y,  z  Is  a  right-handed 
coordinate  system  embedded  in  the  rigid  MUR  beam,  x  coincides  approximately  with 
a  neutral  axis,  x  and  z  are  parallel  to  the  aircraft  symmetry  plane,  X,  Y.  Z  Is 
parallel  to  x,  y,  z  and  embedded  in  the  rigid  aircraft  frame  with  origin  at  the 
initial  rigid,  c.g,  locution  of  the  More  being  released,  x\  y\  //  is  a  right-handed  set 
of  store  body  coordinates. 

Major  assumptions  and  features  of  the  model  are  enumerated  below:  they 
render  the  problem  tractable  and,  hopefully,  physically  representative, 

(1)  Wing,  pylon,  hanger,  shoulder  pad,  and  ejector  unit  deformations  are 
neglected. 

(2)  Coreolls  and  centripetal  accelerations  Induced  by  deformation  are  neglected, 

(3)  “Bernoulli"  beam  pitch  and  yaw  deflections  are  allowed  for  the  MUR 
beam,  Torsional  deformation  is  assumed  to  Induce  no  longitudinal  stresses, 
Longitudinal  oscillations'  arc  neglected. 

(4)  Lug  connections  from  the  parent  rack  to  the  MUR  are  ussumed  to  be 
pinned,  Parent  rack  sway  brace  reactions  are  assumed  to  be  spring  forces. 
It  to  assumed  that  sway  braces  do  not  break  contact  with  the  MUR  during 
deformation  (questionable  for  the  AURO-7A  rock), 
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(5)  Static  and  dynamic  aeroeiastic  forces  on  the  MER  are  neglected  due  to 
lack  of  experimental  data. 

1 6)  Kinematics  of  the  ejecting  foot  and  in-carriage  store  c.g.  are  computed 
from  beam  deflections  and  rotations  at  beam-ejector  unit  connection  cross 
sections. 

(7)  For  a  given  cartridge,  a  single  ejection  force-time  function  is  assumed  for 
all  store  loads.  Ejector  force  dynamics  '.re  assumed  to  be  independent  of 
beam  motion. 

(8)  Beam  motion  acceleration,  aircraft  acceleration,  ejection  force,  and 
aerodynamic  loads  for  in-can iage  stores  are  translated  into  beam  loads  by 
statically  determinant  computations. 

(9)  For  the  application?  given  all  stores  in-camage  are  identical. 

A  short  s-ick  release  of  stores  will  not  deform  the  pylons  ano  wings  greatly. 
For  a  larger  number  of  stores  released  this  assumption  must  be  reexamined.  Some 
of  the  assumptions  involving  structural  deformations,  mass  distributions,  and 
.'.inematics  cculd  he  improved  if  supporting  experimental  data  were  available. 

A.  Influence  Coefficients 

The  MFR  beam  we:ghs  about  210  pounds.  More  than  half  of  the  beam 
weight  is  concentrated  betv,.en  the  ejector  unit  support  cross  sections.  Flexibility  is 
more  important  for  heavier  stores.  Hence,  the  ratio  of  effective  beam  to  store  mass 
is  expected  to  be  1/5  or  less  and  therefore,  the  beam  mass  distribution  need  not  be 
specified  with  great  accuracy.  Therefore,  only  deflections  at  a  relatively  f.-w  cross 
sections  such  as  those  enumerated  1-10  in  Figure  1  need  be  considered. 

For  the  AERO-7A  parent  rack  at  the  sway  brace  cross  sections 

F,  =  -kf  (1) 

For  unit  forces  placed  at  a  collocation  point  (not  including  lugs  where 
deflections  arc  zero)  deflections  are  computed  at  all  collocation  points  (force 
flexibility  coefficients).  For  unit  couples,  placed  at  ejector  unit  and  lug  connection 
cross  sections,  deflections  are  computed  at  the  collocation  points  (moment  flexibility 
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coefficients).  Additional  assumptions  arc  neglect  of  shear  deflection  and  a  constant 
"slue  of  second  moment  of  area  between  collocation  points  or  between  collocation 
points  and  lugs.  Couples  are  due  to  horizontal  forces  and  contribute  less  to 
deflections  than  vertical  forces. 

A  general  equilibrium  equation  may  be  written  as 

o 

p,  and  Mk  include  inertial,  gravity,  maneuver  and  aerodynamic  forces. 

B.  Natural  Frequencies  and  Normal  Modes 

p/  is  now  decomposed  into  distributed  and  concentrated  forces.  The  integral  for 
distributed  forces  is  replaced  by  a  numerical  approximation  (trapezoidal  ntle  is 
adequate ). 


f, 


2  C 
1=1.1 


fij 


-m 


»  dt2 


2 

k  =  1  .6 


CMik  Mk 


(3) 


The  set  of  equations  (3)  are  known  as  collocation  equations  (see  for  instance 
Reference  4).  i  or  j  indices  refer  to  collocation  points.  The  aft  ejector  unit  support 
points  are  at  i  =  j  =  n,jt,  or  k  =  1.2.  Lug  connection  points  are  at  k  =  3.4.  The 
forward  ejector  unit  support  points  are  at  i  =  j  =  n5.nfe  or  k  =  5.6.  ntj  are 
effective  concentrated  masses  at  the  collocation  points.  are  net  vertical  forces  at 
a  cross  section.  Mk  are  net  couples  and  are  due  to  horizontal  forces. 

It  is  convenient  to  choose  an  unloaded  horizontal  beam  subject  to  only  its  own 
weight  as  a  reference  condition.  F  may  be  further  broken  down  into  concentrated 
and  distributed  contributions 


F7j  =  nulgd  -  cos  0Q)  +  0VQ  ]  +  F.^  (4) 

Fzj  is  zero  except  at  the  ejector  unit  support  points.  The  concentrated  forces  at  the 
suppoit  points  arc  due  to  loads  on  the  individual  stores. 
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In  Figure  2(a)  is  shown  the  lorce  system  assumed  to  act  on  a  centerline  store 
and  hanger  supports  removed  as  a  free  bod>  from  the  MER  beam.  The  hangers  and 
ejector  unit  are  not  shown,  only  the  forces  and  points  of  ac'ior 

The  assumption  of  a  pinned  connection  with  equal  horizontal  forces  at  the 
hanger-beam  connection  point  ienders  the  computation  of  f^(.  f  and  f 
siaticalh  determinant.  From  moment  and  force  eoailibrium 


f  =  r  /n 
\  c  X  c '  * 


',u  =  <FZcXf  +  MYc+FXcZc>/Xh 
'/vl  =  (F7.cXa  "  ^Yc'  FXcZc^Xh 


(5) 


In  Figure  2(b)  is  shown  a  sideview  of  the  force  system  acting  on  a  shoulder 
station  store.  Again  from  equilibrium 


*;<  -  Fxs/4 


7  Sl* 


|Fzsxr  +  FxA“  Az/2)  +  MYll/xh 
|FZsxa  -  FXs(zs  -  Az/2)-  My  J/xh 


(6) 


Support  point  net  forces  and  moments  are  given  by  equations  like 


F  =  f  +2f 

L  a  7C  3  7Sa 


=  Mf  *  z,fx,. +(2z2+Az)fxs 


(7) 


After  ejection  is  initiated  the  centerline  store  forces  of  a  separating  store  are 
replaced  by  the  ejection  force  and  (5)  is  replaced  by 


f  =  0 


-Fe(t)(xf-  xc  )/xh 

-Fc(tKx3  +  V/xh 
6 


(8) 


. . .  am  mtfiu^lliMiililll^ 


The  couples  Jt  the  Lius  arc  given  by 


M.  =  M  .  = 


'.I^XI  +  2fpx  3  +  FX4>+FX2* 


(9) 


Thus  far  the  sets  ot  conations  (6)-(9)  are  valid  for  all  store  loads.  Some  of  the 
stores  and  hence  loads  ma>  he  absent. 

Inertial  loads  at  the  store,  c.g.  must  be  defined  in  order  to  complete 
formulation  of  the  characteristic  equations. 

It  is  assumed  that  the  store  and  ejector  unit  assume  the  angular  attitude  of  the 
secant  of  the  deflection  curve  between  support  cross  sections.  For  a  centerline  stole 


"  'S  +  V 


Z  e 


-  Nntxjf +  -\ia> 

(Xa  +  Xf) 


M 


-  ~  B,£.  -  ft* 


Yc 


(.V3  +Xf) 


(10) 


Shoulder  station  equations  are  similar  with  z  replaced  b>  z^.  The  horizontal 
acceleration  contributes  least  to  the  natural  frequency. 

Substitution  of  the  sets  of  equations  (6KI0)  into  equation  (3)  with  only 
inertial  forces  acting  yields  a  homogeneous  set  of  equations  of  the  form 


f, 


V 

J=l  A 


Knm, 


dt2 


1 1 ; ) 


Tlie  eigenvalue  problem  is  given  by 


XZ.  =  f 


112) 
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k,  i»  no!  symmctru  and  heme  tin-  tor.nulation  is  non  self-adjoint  The 

vflJlJCJcnMK  CijUil  lion  coefficients  are  obtained  by  application  of  the 

(V-lcy-Hamlton  theorem  i  Reference  5i  A  good  first  iteration  for  the  dominant 
ciccnv.iiuc  i>  siv.n  In  Sv Ivcster's  theorem  (Reference  5).  The  eigenvector  for  the 


r.-nvuhie  X  is  given  b\ 

n 


X  /.  -  ZK  .,iZ 

r  in  it  •  jn 


(13) 


liu  adjoint  eigenvector.  /  lor  X„  is  given  bv 


x  z  -  >;k  in  z 

n  in  ||  j  in 


(  14) 


Orthogonality  i*  gnen  bv 


imZ  Z„  =  0 

i  in  (in 


(15) 


C.  Dynamic  Response  of  MKR  Beam  to  Stou  Release 

The  reference  condition  from  which  the  dellection  is  computed  is  the 

horizontal  Ml-R  beam  under  its  own  weight.  Normally  only  in-carriage  airloads  are 
available  (airborne  or  wind  tunnel)  for  the  store  about  to  be  released.  However, 

ill-carriage  airioads  for  all  stores  for  a  given  release  condition  were  available  for  the 

A-?!)  from  LTV  for  the  M-117  bomb  (Reference  6).  If  only  airloads  are  available 
lor  the  store  about  to  be  released,  the  reference  condition  must  be  the  deflected 
position  just  prior  to  release. 

Incremental  loads  due  to  static  aeroelastic  deformation  are  not  considered. 
However,  the  deformation  is  considered  when  the  initial  angle  of  attack  of  the 
released  store  with  the  free  stream  is  computed  for  the  ejection  phase.  It  is  assumed 
that  during  the  ejection  phase  the  airloads  on  stores  remaining  in-carriage  do  not 
change.  Static  measurements  of  in-carnage  airloads  indicate  fairly  large  dif**:re»ces 
between  a  stvbomb  and  five-bomb  configuration  or  more  generally  between  an 

VHomb  and  N-l-homb  configuration.  Conceivably,  transient  in -carriage  loads  could  be 
measured  by  in-carriage  balances.  These  dynamic  loads  would  be  due  to  MLR  beam 
vibrations  and  released  store  tnuiMti  ’n  through  the  imciference  near-flow  field. 
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Without  any  cxpen.nenui  data  available  the  transient  between  static  conditions 
would  have  to  he  assumed  Such  considerations  have  r.ot  been  applied  h  re. 


The  general  equations  of  motion  may  be  w'itten  as 


•A.K'.m.5c  •t"F."hL. 

Lt  is  the  -'"latiop  due  to  cons*ant  air.  gravity,  and  maneuver  loads. 
\m  eiy  iivector  expansion  is  assumed  (truncated  at  N  n  odes) 


tlb) 


n  -  1  .N 


(lit 


Substitution  of  (17)  into  (16)  with  application  of  (15)  yields 


T  +  cc;T_ 

rs  n  ~ 


M. 


*  •»1Z1„L,-Fe(«>  *  "AA. 

i ,i  i- 1  ,t 


(IS) 


1 18)  is  readily  solved  as 


T„(0) 


_  „  .  mtZin L.(  1  -  cos t ) 

T  =  T  (O)cosco  t+  -  sin  u;  t  +  Z  - - 

«„  "  «»!.*  M„ 


i-7  v  m«KA"  f  S'n  Wn(t  '  T)Fc(ThJr 


(19) 


f.lOirrijZ, 


TjO)  =  Z 
n  i  .v  M 


f.(0)m.Z._ 

r  (0)  =  z  - — — 

n  Ni 


CO) 


iqu.it ions  si 9)  and  (17)  give  Hie  beam  response  to  store  release  and  ejection  force. 
it<Ui  Is  non  zero  only  tor  the  ripple  case.  Fhe  response  is  computed  for  about  100 
mi!ii>cvuiid'.  and  stored  tor  later  combination  with  store  dynamics  during  the 
election  phase 


HI.  EJECTION  PHASE  STORE  DYNAMICS 


A.  Ejection  Kinematics 

Kinematics  of  the  store  and  »>Sior  unit  ar-.  computed  relative  to  th*'  X.  Y.  7 
coordinate  s\  .tem  \c-.isun  eei-met r\  is  shov*r.  in  E'igure  3. 

f  and  a  arc  the  same  points  as  shov. n  in  Figure  2.  0  and  0  aie  small  angles, 
p  is  a  point  on  the  bottom  sin  face  of  the  ejector  unit  where  the  piston  foil  exits. 
p_  .s  ’he  rigid  case  location  of  point  p.  pp'  remains  perpendicular  to  the  bottom 
surface  o!  the  cje.tor  unit 

The  point  p  coordinates  are  given  as 


X 

p 


x  +  / 

c  c 


=  X  +  l  O 

V  C 


z 


p 


(2!) 


Initial  conditions  for  the  store  are  given  b> 


0  (0)  =  0(0) 

\ 

X%(0)  =  Os 

Z  ;0)  =  f  (0)  -  (x  -  \  kxO) 

0  (0)  =  0(0) 

> 

U  <0)  =  z  o<0) 

s  c 

W  (0)  =  t  (0)-  (X  -  x  >ctf0)  <22 > 

'  ^3  p  v 


Subsequent  values  of  X4  and  7.s  arc  governed  by  the  dynamics  of  store  separation. 


The  length  ot  the  piston  throw  tpp')  i>  given  by 


L) 

7  -  7  +  0  (X  -  X  > 

s'*  p  VS  p 


The  lor«.e  moment  arm  is  given  b> 


-  v  ,:4‘ 

B.  lijection  Fquations  of  Motion  :n  X.  Z  Coordinates 

In  figure  4  is  shown  some  additional  velocity  vector  and  geometry  information. 
The  equations  o!  motion  in  body  coordinates  are 


M  tu'  +  q’w'i  =s  -pV*  SC.  -  M  e  sin  d  +  F  ttxo-0  ) 

%  1  ■>  r  A  c  » 


M  iw’  -  q’u'i  =  -  Q'SC  s  +  M  g  cos  0  +  Ffit ) 


Bq'  =  O’SD 


C  * 

1  m  '  2V'  ^  ci»! 


“  v  F  (t) 

tr,  c 


(25  > 


During  ejection  the  aircraft  velocity  magnitude  remains  at  VQ.  The  aircraft  is  pulling 
up  at  a  constant  angular  rate  6  It  is  assumed  th3t  the  angle  of  attack.  a0.  of  the 
X  axis  remains  constant 

Velocities  of  the  store  relative  to  the  X.Z  coordinates  3re  introduced 


=  Vocosao+07-o  +  Z,0  +  Uv 


=  V0  sin  0^ 


0xn  -  X  6  +  W. 


(26) 


From  kinematics  one  can  obtain  in  addition 
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^  AIRCRAFT  REFERENCE  LINE 


Klnwruttlc  G#om»try 


p  *  p  0 


V'2  =  u'2  +  W'2 


V2 

*  A 


1  +  2  cos 


a„ 


(Oz  o  +  U4) 


V„ 


V2[]+0«.01)]  -  v2 


V' 

Ma  =  —  «*  Man 


O’-Oo 

w'  Wc  -  0xo 

a- ocQ  -an*0+  — 5 - - 

0  U  0  s  y 

v0 

sin  0  ^  sin  0Q  +  cos  0Q(0t  +  0^) 
cos  0  =s  .-os  0Q  -  sin  0p(0t  +  0 s) 


(29) 


In  addition 


ax  =  W  0  +  Us  +  0w 

az  =  +Ws  -  0u  (3Q) 


liquations  (26H30)  arc  applied  to  the  set  of  equations  P5)  with 
order  terms. 


retention  of  first 


0,  -  -«2W,  ♦  V0a0  MEM  .  6 sin  „  M  .  OoSCn. 

0  M.  M 


W  =  0V  +  ft  -  2ol 
°  Ms  Ms 


r  dCS  3Cn 

+  —  Zs 


+  g  cos  0. 


BO  = 


Q0SD 


z  vo 


-FclXp-Xi  +  «Z,-Zl| 


p' 

(31) 
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Solution  of  equation  (31)  begins  with  initial  conditions  (22).  Computation  ends 
when  the  sector  has  reached  its  end  of  stroke. 


fe,  =  £ 

pt  c 


(32) 


The  ejection  velocity  is  defined  as  the  store  c.g.  velocity  vector  at  the  end  of 
stroke  relative  to  its  initial  velocity  vector.  Components  of  the  ejection  velocity 
along  and  perpendicular  to  the  initial  velocity  vector  are 


Wc  *  ws  -  0Vote 

uc  -  us 


Wf  is  normally  what  is  called  the  ejection  velocity  since  Ue  <  VQ. 


(33) 
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IV.  APPLIED  COMPUTATIONS 


The  model  developed  in  the  sets  of  equations  ( 1 )  - (33)  was  programmed  for 
the  CDC  6700  with  general  input  variables.  Specific  applications  are  indicated  below. 

A.  MLR  Structural  Dynamic  Data 

Structural  and  geometric  c.jta  were  taken  from  References  7-11.  The 
configuration  chosen  was  a  MER  rack  shifted  completely  forward.  Parent  rack  was 
the  MAU-9  (lugs  a.td  sway  braces  essentially  in  the  same  plane). 

Light  collocation  points  are  taker,  at  x  locations  x  =  0.,  1.08,  2.77,  4.65. 
5.48.  7.60.  8.90.  and  10.58  feet.  Ejector  unit  connection  points  are  at  x  =  1.08, 
2.77,  8.90.  and  10.58  feet.  Lug  connections  are  at  x  =  3.81  and  6.31  feet.  For 
collocation  points  on  the  strongback,  I  was  computed  as  43.1  in4:  for  other  points 
I  =  27.7  in4.  The  elastic  modulus  for  aluminum  alloy  7075  is  E  =  10.4  X  106  psi. 
Effective  weights  at  collocation  points  (in  corresponding  order  to  x  points  above) 
were  estimated  as  5.43,  36-01,  36.79,  9.94,  9.94,  15.41,  38.28  and  34.13  lbs.  Other 
pertinent  geometric  data  are  Uc  =  .35.  ze  =  .479,  z£  =  -.259,  zf  =  .252, 
z,  =  -.198.  Az  =  .2375.  xp  =  .6875.  xo  =  5.,  and  xh  =  1.6875  ft.  The  remainder 
of  the  data  are  dependent  on  the  configuration. 

One  check  on  the  structural  model  is  a  comparison  of  a  nose  deflection 
computation  with  data  from  Reference  11.  The  configuration  is  3  M-117  bombs 
forward  and  2  shoulder  M-117  bombs  aft. 

Test  loads  are  taken  from  Reference  10.  Loads  on  the  forward  stores  are 
F 2 =  Fz,  =  6.57  X  103  lbs,  FXc  =  FXl  =  -1.76  X  103  lbs.,  and 

MVc  =  Mys=  5.31  X  103  in-lbs.  Loads  on  aft  stores  are  FZs  =  6.54  X  103  lbs. 
FXs  =  -1.76  X  103  lbs.  and  Mys  -  7.61  X  103  in-lbs.  Geometric  data  from 
References  7  and  12  are  xc  =  .025,  zf  =  1.25,  zs  -  .703,  and  D  =  1.333  ft. 
Computation  yields  2.8  inches  for  the  nose  deflection  as  compared  tc  3.1  inches 
from  Reference  11.  This  indicates  that  the  structural  model  is  at  least  representative 
of  the  MLR  in  pitch  plane  bending. 

The  ejector  forcing  function  used  is  a  mean  curve  for  data  for  the  Mk  2-1 
cartridge  taken  from  Reference  8.  The  forcing  function  compared  with  the  original 
data  is  reproduced  in  Figure  5. 
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B.  Applications 


Sample  computations  were  made  corresponding  to  0Q  =  -30°.  \1a  =  .9, 
and  an  altitude  of  10000  ft.  Corresponding  air  densitites  and  speed  of  sound  are 

based  upon  the  ICAO  standard  atmosphere.  The  pull-up  rate  was  varied  between 

0  =  0  and  0  =  .1  rad/sec  (about  a  3.9  g  normal  acceleration). 

As  indicated  previously,  complete  in-carriage  aerodynamic  data  were 
available  only  for  the  M-117  (103  tail)  bomb.  Store  inertia!  data  from  Reference  12 

are  Ms  =  824  lbs.  and  B  =  1609  lb-ft2.  Aircraft  aerodynamic  data  were  for  an 

A7-D  aircraft.  The  rack  was  located  on  the  center  wing  pylon.  Neighboring  pylons 

were  unloaded.  Inboa.J  and  outboard  shoulder  in-carriage  pitch  plane  aerodynamic 

loads  were  averaged  to  obtain  equivalent  pitch  plane  loads.  The  initial  angle  of 
attack  of  the  aii  -raft  reference  line  was  assumed  to  vary’  linearly  with  0  from  3.5° 

at  0  =  0  to  10°  at  0  =  .1  rad/sec.  Angles  of  attack  of  the  store  are  obtained  by 

subtracting  3°  (see  ftr  instance  Reference  13). 

Beam  response  was  computed  for  a  6-bomb  configuration  and  an  aft 

centerline  release  with  0  =  .1.  Because  of  the  short  distance  from  the  rear  lug  to 
the  ejedor  unit,  deflections  were  small  and  induced  velocities  at  the  forward 

centerline  store  c.g.  were  small  (.2  ft/sec  maximum).  Hence,  flexibility  for  aft 

cenit-rlinc  releases  was  considered  negligible.  Induced  ripple  effects  were  also 

condudet  to  be  negligible  for  this  case.  The  remaining  computations  were  made  for 
a  5  bomt  configuration  a.id  a  forward  centerline  release. 

Free  field  aerodynamics  were  obtained  from  References  12.  14  and  15. 
Cn-q  and  CmQ  were  estimated  from  local  slopes  at  the  initial  angles  of  attack. 

Three  modes  were  U3vj  to  compute  the  MER  beam  response. 

For  the  M-l  1 7  bomb  computations  were  made  for  3  cases 

(1)  Flexible  MER  with  variable  0. 

(2)  Rigid  MER  with  variable  0. 

(3)  Rigid  MER  with  0  =  0  but  the  came  ••  .ual  angles  of  attack  as  for  (1) 

and  (2). 

Nominal  ejection  velocity  for  the  M-117  from  Reference  8  is  about  6.4  ft/sec. 
First  and  second  mode  frequencies  were  computed  as  44.3  and  68.9  rad/ see.  In 
Figure  6  is  plotted  the  flexible  MER  ejection  velocity,  Wc  (Case  (1)).  and 
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(  FT / SEC ) 


contributions  due  to  maneuver  and  flexibility.  Contribution  due  to  pull-up  rate  is 
labeled  (AW  )  =  W  ,  -  W  Contribution  due  to  flexibility  is  labeled 

iAW),  =  Wc,  -  \V  Figure  7  is  a  plot  of  angular  rate.  qc  (case  (1)). 

(Aq', >,;  =  ,  -  q'3  and  (Aq')f  =  o',  -  q',.  The  numerical  subscripts  refer  to  the 

case  numbers  above . 

For  the  M- 117  flexibility  decreases  the  ejection  velocity  due  to  aerodynamic, 
ejector,  and  gravity  forces  by  some  8.7  percent  at  0  =  0  and  lesser  amounts  for 
increasing  pull-up  rates.  Note  that  the  effect  of  pull-up  rate  alone  is  significant.  For 

abom  .017  rad/sec  maneuver  and  flexibility  effects  are  about  equal.  Figure  7  shows 

that  the  ejection  angular  rate  of  the  store  is  affected  to  a  much  greater  extent  by- 
flexibility  (30  percent  for  Q  =  .1).  In  this  case,  the  angular  rate  magnitude  is 
decreased  or  the  effect  of  initial  disturbances  is  alleviated. 

Computations  were  also  made  on  a  hypothetical  300  lb.  1  omb  in  order  to 

investigate  effects  of  flexibility  and  maneuver  on  a  more  lightly  loaded  beam.  The 
bomb  was  assumed  to  have  the  same  aerodynamic  coefficients  and  values  of  x  as 
the  M- 117  bomb.  Input  variables  which  are  different  than  the  M- 117  case  are 
B  =  600  Ib-lt2.  D  -  .75  ft.  /  =  .058  ft.  and  /  =  .497  ft.  From  Reference  8  {he 

i.  > 

nominal  ejection  velocity  is  about  9.2  ft/sec.  Frequencies  of  the  first  two  modes  are 
71.9  and  113  rad, sec. 

Figures  (8)  and  <9i  show  that  the  maneuver  contribution  is  still  significant  for 
the  more  lightly  loaded  MhR  beam. 
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300  Lb.  Store  Ejection  Velocity 
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V.  CONCLUDING  REMARKS 


A  model  has  been  developed  which  gives  reasonable  estimates  of  the  effect  of 
maneuver  and  MER  flexibility  on  store  ejection  conditions  for  pitch  plane  motion. 
It  appears  that  the  ejection  velocity  is  decreased  by  no  more  than  iO  percent  for 
the  heaviest  stores.  However,  the  ejection  angular  rate  due  to  interference 
aerodv namics.  ejector  force,  and  gravity  is  significantly  affected. 

In  future  work,  it  is  hoped  that  results  will  be  substantiated  by  experimental 
data  on  store  k. nematics  ar.d  MER  structural  dynamics. 

It  is  planned  to  extend  this  stud}  to  yaw  plane  and  torsional  deformation  and 
Mx-deeree-of-freedom  motion.  Yaw  plane  values  of  stiffness  are  about  half  that  of 
the  pitch  plane.  Parent  rack  swa\  bracing  is  relatively  weak  (especially  for  the 
Al-RO-7  rack).  Hence,  it  is  anticipated  that  the  effect  of  MER  deformation  is  more 
important  for  shoulder  station  releases. 
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APPENDIX  A 
LIST  OF  SYMBOLS 


fv 


LIST  OF  SYMBOLS 


B  Pilch  moment  of  inertia  of  store  (slug-ft2) 

c  Speed  of  sound  (ft /sec) 

CA  Axial  force  coefficient 

CD  Drag  coefficient 

CD0  In-carriage  drag  coefficient 

Cf  Force  fiexibility  influence  function  (ft/lb) 

Cfjj  Force  fiexibility  coefficient  matrix  (ft/lb) 

C  Pitch  moment  coefficient 

m 

C  „  In-carriage  value  of  C 

mi»  m 

C  .  Pitch  damping  coefficient 

CM  ik  Moment  flexibility  coefficient  matrix  (lb~ 1 ) 

CN  Normal  force  coefficient 

CN0  In-carriage  value  of  CN 

D  Store  reference  diameter  (ft) 

Fc  Ejection  force  (lb) 

f  x  component  of  force  reaction  at  hanger-beam  connection  point  for  a 

centerline  MER  station  (lb) 

Fx  X  component  of  force  acting  at  centerline  store  c.g.  (lb) 
f  Similar  to  fxc  for  shoulder  station  (lb) 

Fv  Similar  to  Fv„  for  shoulder  station  (lb) 

Fx,  FXc  for  aft  station  (lb) 


A-l 


F\  2  FXc  for  forward  station  (lb) 

Fx3  FXs  for  aft  station  (lb) 

I'x  4  FXs  for  forward  station  (lb) 

Fz  Concentrated  part  of  Fz  (lb) 

Fz  z  component  of  total  force  acting  at  a  MER  beam  cross  section  (lb) 

Fza  Fz  at  an  aft  end  of  an  ejector  unit  (lb) 

FZc  Z  component  corresponding  to  FXc  (lb) 

f  z  component  of  force  corresponding  to  f  for  aft  connection  (lb) 

f  ,  Similar  to  f  for  forward  connection  (lb) 

FZj  Fz  at  jth  collocation  node  (lb) 

F  .  T,  at  jth  collocation  node  (lb) 

FZs  Similar  to  FZc  for  shoulder  station  (lb) 

fZ4a  Similar  to  fzca  for  shoulder  station  (lb) 

fzsf  Similar  to  fzcf  for  shoulder  station  (lb) 

g  Acceleration  of  gravity  (ft/sec2) 

I  Second  moment  of  area  of  MER  beam  cross  section  (ft4) 

k  Sway  brace  spring  constant  (Ib/ft) 

Kci  Flexibility  matrix  associated  with  ejector  force  (ft/lb) 

Ky  Flexibility  matrix  associated  with  beam  acceleration  (ft/lb) 

C  Beam  length  (ft);  number  of  beam  collocation  points 

8  ,  Piston  throw  length  (ft) 
pi 


A-2 


Li 

Ma 

Ma() 


M 


a 


M, 


m 


M 


k 

Mv, 

M*r 

M„ 

Ms 

Mv  . 


M 


Ys 


N 


P, 

q' 

0 

Oo 

s 


Deflection  vector  associated  with  time  independent  forces  (ft) 

Mach  number 

Initial  Mach  number 

Couple  at  aft  ejector  support  (ft -lb) 

Couple  at  forward  ejector  support  (ft-lb) 

Effective  mass  at  ith  collocation  node  (slugs) 

Couple  at  kth  node  (ft-lb) 

Similar  to  Mk  for  aft  lug  cross  section  (ft-lb) 

Similar  to  Mta  for  forward  lug  cross  section  (ft-lb) 

\-  I 

Store  mass  (slugs) 

Moment  acting  on  centerline  store  (ft-lb) 

Similar  to  My  for  shoulder  store  (ft-lb) 

Number  of  modes  in  J.  expansion 
Load  distribution  on  MER  beam  (lb/ft) 

Angular  rate  of  store  (rad/sec) 

Dynamic  pressure  (lb/ft2) 

Initial  0  (lb/ft2 ) 


t  Time  (secs) 

t  Ejection  time  (sees) 
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Vo 

w 
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W 


W 

c 

Ws 

x 

x' 

X 


X. 


Time  function  associated  with  nth  mode 
X  component  of  velocity  of  store  (ft/sec) 
x'  component  of  store  velocity  (ft/sec) 

Component  of  ejection  velocity  parallel  to  initial  aircraft  velocity  (ft/sec) 

X  component  of  velocity  of  store  relative  to  X,  Y,  Z  origin  (ft/sec) 

Store  velocity  magnitude  (ft/sec) 

Aircraft  velocity  magnitude  (ft/sec) 

Z  component  of  velocity  of  store  (ft/sec) 
z'  component  of  store  velocity  (ft/sec) 

Similar  to  Uc  for  perpendicular  component  (ft/sec) 

Similar  to  Us  for  Z  direction  (ft/sec) 

Longitudinal  axis  along  MER  (ft) 

Longitudinal  store  body  axis  (ft) 

Axis  parallel  to  x  with  origin  at  initial  rigid  c.g.  position  of  store  (ft) 
Initial  x  distance  from  store  c.g.  to  aft  ejector  unit  support  (ft) 

Similar  to  xa  for  forward  support  (ft) 

Initial  x  distance  from  piston  to  rigid  store  c.g.  (ft) 

=  xa  +  xf 

x  distance  from  aft  support  to  ejector  foot  (ft) 

X  location  of  point  on  bottom  surface  of  ejector  unit  adjacent  to  piston 
foot  (ft) 
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X  location  of  store  c.g.  (ft) 

x„  x  distance  from  aircraft  c.g.  to  in-carriage  rigid  c.g.  of  store  (ft) 

z  Axis  perpendicular  to  x  (ft) 

Axis  perpendicular  to  x'  (ft) 

Z  Axis  perpendicular  to  X  (ft) 

z  z  distance  from  bottom  of  MER  beam  to  initial  store  c.g.  (ft) 

zv  z  distance  from  bottom  of  MER  beam  to  bottom  of  ejector  unit  (ft) 

nth  mode  vector 
nth  adjoint  mode  vector 
/.v  7  distance  from  neutral  axis  to  lug  (ft) 

Zp  Z  component  corresponding  to  Xp  (ft) 

z()  z  correspondent  to  xQ  (ft) 

7 j  Distance  from  neutral  axis  to  bottom  of  MER  beam  at  ejector  unit  (ft) 

7 2  Similar  to  z( .  but  distance  to  upper  shoulder  support  (ft) 

a  Store  angle  of  attack  (rads) 

a0  Initial  value  of  a  (rads) 

Az  z  distance  between  shoulder  supports  (ft) 

f  Beam  deflection  at  aft  support  (ft) 

Beam  deflection  at  forward  support  (ft) 

Deflection  at  ith  collocation  node  (ft) 

0  Store  elevation  angle  above  horizontal  (rad) 
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PQ 
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Constant  pull-up  rate  of  aircraft  (rad/scc) 
Initial  elevation  angle  of  store  (rad) 

Angle  between  x'  and  X  axis 
Eigenvalue  (sec/rad) 
nth  eigenvalue  (sec/rad) 

Air  density  (slug/ft3) 

Initial  p  (slug/ft3) 

-  -  (rad) 

xh 

Frequency  of  nth  mode  (rad/sec) 
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